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Abstract. We study the spectral properties of the transfer matrix for a gonihedric 
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1 Introduction 



The triviality of simple random surface models on the lattice, due to the dominance 
of branched polymers |l|], and the non-scaling of the string tension in analogous 
Euclidean invariant models based on dynamical triangulations made it imperative 
to study models with an action depending on the extrinsic curvature of surfaces. For 
a review of work in this area up to 1997, see 0. 

One of the models with extrinsic curvature dependent action is the so-called 
gonihedric random surface model introduced in [|J and further studied in [|, ||, 
[7j. This model is unstable in the simplest cases ||. However, there is a lattice 
discretization of the model || |10[ which is not plagued by stability problems and is 
more amenable to analytical study, see |TTJ-|jl6[. Another lattice model with similar 



properties was introduced and studied in |fL7|. 



In |18[ a simplified version of the three-dimensional gonihedric lattice model was 
studied and the two largest eigenvalues of the transfer matrix were determined. In 
this paper we finish the calculation of the spectrum of the transfer matrix for the 
model and find all the eigenf unctions explicitly. 



2 Transfer matrix for loops 

Let T3 denote a sublattice of Z 3 of size MxMxN with periodic boundary conditions. 
We shall think of the third coordinate direction where the lattice has extension N 
as the 'time' direction. The configuration space of the system we wish to study is 
a collection of subsets of the plaquettes in T 3 , denoted At, which we refer to as 
singular surfaces or simply as surfaces. A collection of plaquettes M belongs to Ai 
if and only if any link in a plaquette in M belongs to an even number of plaquettes 
in M. This means that the surfaces in At do not overlap themselves but they can 
intersect themselves at right angles and they are closed. Note that the surfaces need 
not be connected. It is not hard to see that the conguration space Ai is identical to 
the configuration space of the three-dimensional Ising model on the periodic lattice 
T 3 . The surfaces are in one to one correspondence with phase boundaries. 

We say that a link £ in a surface M e Ai. is an edge if i is contained in exactly 
two plaquettes that meet at a right angle. We denote the collection of all edges in 
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M by S(M). The action S(M) of the surface M is defined as the total number of 
edges in M 

S(M) = #£(M). (1) 

Note that no action is associated with those links where the surfaces crosses itself. 
The partition function is given by 

Z{f3) = e~^ M \ (2) 
MeM 

We define an intermediate plane in R 3 to be a plane that lies parallel to but in 
between planes where one of the coordinates takes an integer value. If we slice 
a surface M (regarded as a subset of R 3 ) by an intermediate plane we obtain a 
collection of links which we shall call generalized loops or simply loops. The collection 
of all possible loops that can arise in this way will be denoted n. A collection of 
links P in an intermediate plane belongs to II exactly when each vertex in P belongs 
to an even number of links in P. We note that the loops II coincide with the phase 
boundaries that arise in the two-dimensional Ising model. It will not come as a 
surprise that this close connection to the Ising model will make the present model 
exactly soluble in a very strong sense. 

For P G II we let \P\ denote the number of links in P. Note that if Pi, P 2 G IT, 
then the symmetric difference 

(Px u p 2 ) \ (Pi n p 2 ) = p ap 2 g n (3) 

is again a loop in IT. It is easy to check that 

|PAP 2 | = | Pi | + |P 2 | -2|pnP 2 | (4) 

and d(Pi, P 2 ) = |Pi AP 2 | is a metric on II. 

Let k(P) denote the number of corners of P, i.e. the number of vertices of order 2 
whose adjacent links are at right angles. Let us denote by Ti the space of real- valued 
functions on II. Let Kp be the operator on TC whose matrix (kernel) is given by |TT 



K P {P U P 2 ) = exp + k{P 2 ) + 2|PiAP 2 |)) . (5) 

With this notation the partition function of the model can be expressed as 

Z N {(3) = Tr Kg (6) 



by slicing surfaces in intermediate planes orthogonal to the time direction. This 
partition function can be evaluated if we ignore the interesection term \P X n P%\ 



in the exponent of Kp or if we replace k(P\) + k(P2) by k(PiAP2) [IT]. In this 
approximation the model becomes a stack of noninteracting two-dimensional models. 
The critical beahaviour of the model is identical to that of the two-dimensional Ising 
model and the free energy can be computed exactly by standard methods |19| . 



Here we will make an approximation different from the one of [11 1. We drop the 



curvature terms in the action and study the simplified transfer matrix 

KpiPi, P 2 ) = exp {-2/3(\P 1 AP 2 \ - M 2 )) , (7) 

where we have inserted a normalization factor M 2 in the exponent for later conve- 
nience. Dropping the curvature term is the same as neglecting the action associated 
with edges in the time-direction. We note the mathematical analogy of the transfer 
matrix (^) with the transition function exp(— j3\x — y\), x,y G R d , for random walk 
in R d . This makes it reasonable to regard Kp as describing the diffusion of loops. 

With abuse of notation we denote the new transfer matrix (0) by the same 
symbol as the old one ([5]). This should not cause any confusion since we stick to the 
notation of Eq. (|7|) in the sequel. 

3 Eigenvalues of the transfer matrix 

In this section we solve the eigenvalue problem 

£ k p (p, Q)UQ) = HP)Mn (8) 

Q&n 

ipi € TC. We prove that only the eigenvalues depend on f3, not the eigenfunctions. 
The eigenvalues can all be expressed in terms of the partition function, the internal 
energy and the spin correlation functions of the two-dimensional Ising model. The 
eigenfunctions are explicit functions on II which can be normalized to take only the 
values 1 and —1. 



Note first that 



Q Q 

= Ao(/9), (9) 
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where we have in the first step shifted the summation variable from Q to PAQ 
(permissible since the mapping P i— ► PAQ is bijective) and A (/3) is the partition 
function of a two-dimensional Ising model on a periodic M x M lattice. This lattice 
will be denoted T 2 . The loops Q are the phase boundaries of the Ising model and 
the Ising spin variables sit on the lattice dual to T 2 . 

It follows from Eq. (||) that the constant function is an eigenfunction of Kp and 
Aq(/3) is the corresponding eigenvalue. Since all matrix elements of Kp are positive 
and all entries in the eigenvector corresponding to Ao have the same sign, it follows 
from the Perron-Frobenius theorem that A is simple and it is the largest eigenvalue 



of Kp. We conclude |18[ that the free energy per site in the present model is the 
same as the free energy per site in the two-dimensional Ising model. In particular, 
the two models have the same critical point and the same specific heat. However, the 
correlations are different, and we proceed to the calculation of the next eigenvalue 
of Hp. 

Let us introduce the notation 

(P\Q) = 2M 2 -2\PAQ\, (10) 

so Kp(P,Q) = exp (/3(P\Q)). We find it convenient to think of (P\Q) as the 'inner 
product' between the loops P and Q even though the set of loops is of course not a 
linear space. We shall refer to {P\Q) as the invariant product on IT. 

If P is a loop we let P denote its complement, i.e. the loop made up of exactly 
those links in T 2 that P does not contain. We let denote the empty loop and U 
the loop that contains all links in the lattice. Then = U and P = PAU. The 
invariant product is clearly symmetric in its two arguments and it has the following 
properties: 

(P\P) = 2M 2 (11) 
(P|0) = 2M 2 -2\P\ (12) 
(P\Q) = ~(P\Q). (13) 

One can easily construct a finite dimensional inner product space which contains all 
loops in a natural fashion such that the inner product of loops coincides with the 
invariant product and the loops lie on a sphere of radius \[2M. 
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The invariant product is invariant under a nonabelian group Q of motions in 
loop-space generated by 

• the translations of loops 

T a :P^P + a, (14) 

a e T 2 , 

• 'generalized antipodal maps' 

A Q :P^ PAQ, 

gen, 

• reflections and rotations in T 2 . 

Invariance means that (gP\gQ) = (P\Q) for any g e Q. As explained above we can 
view the set of loops as a sphere and we are free to regard the empty loop as the 
north-pole and the complete loop U at the south-pole. Loops increase in size as 
one moves from north to south. The mapping P i— > P is in this picture the usual 
antipodal map and the mappings Aq all have the property Aq = I where / is the 
identity map. 

The group Q acts on the functions in Ti in a natural way: 

g^(P) = ^(g- 1 P). (16) 

Let O n denote the linear operator on H with matrix elements 

O n (P,Q) = (P\Q) n . (17) 

Then O n is invariant under the action of Q on Ti, i.e. 

gO n g- 1 = O n (18) 

and the operators O n commute with each other. Since they are symmetric it follows 
that they have common eigenvectors and the same applies of course to the transfer 
matrix 

oo on 

K ?=Y.ZjO n . (19) 
n=0 n - 

This proves that the eigenvectors of Kp are independent of (3. 
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(15) 



In view of the invariance of the transfer matrix under Q it is natural to look for 
eigenvectors which depend on P only via the invariant product. Let us define 

^q(P) = (P\Q). (20) 

We claim that 

Kfrta = Ai (21) 

for any Q E II, where 

Ai(/3) = ^E^ (P|0> ™, (22) 



2M 2 j, 



so the eigenvalue ratio AiA 1 is minus the internal energy per link of the Ising model 



on T 2 18 



In order to prove Eqs. fl2~T|) and (22) it is convenient to introduce a function rjg 



on II, where £ is a link in the lattice T 2 , defined as 

m(P) = {-\ (23) 

Note that 

(P|Q) = ]>>(PAQ), (24) 

e 

where the sum runs over all links in the lattice T 2 . Consider the function 

: P i-> ru{PAQ) (25) 

on n. lf££Q then £ E PAQ if and only if £ E P. Similarly, if £ E Q then £ E PAQ 
if and only if £ ^ P. We conclude that 

m (PAQ)= Ve (P) m (Q). (26) 

It follows that 

J2e^ V e(PAQ) = £ e^rj e (P) Ve (Q) 
p p 

1 



2M 2 Jt 



J2e mP) (0\P) m (Q). (27) 



The last equality is obtained by observing that the sum 



J2e mP) m(P) (28) 
p 
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is independent of £ due to the translational and rotational invariance of the invariant 
product and using Eq. This proves Eqs. (|2l|) and fl22f) since it suffices to verify 

Y J Kp{P\P)(P\Q) = A 1 (P'\Q) (29) 
p 

for P' = due to the invariance of Kp under Q. In fact we have shown that the 
functions Qq/ are all eigenfunctions of Kp with the eigenvalue Ai. 

We now use the elementary functions $q j to construct the higher eigenfunctions 
of Kp. Let t\ and £ 2 be two distinct links. Define 

®Q,QM{P) = ^ QlA (P)^Q 2 ,eAP)- (30) 

Then, by Eq. @, 

(Kp$ QlQM ) (0) = A 2 (£ u £ 2 )$ QlQ2Ae2 (0), (31) 
i.e. ^QiQa/^a is an eigenfunction of Kp with eigenvalue 

A2(£u£2) = Y,e mP) VtAP)Vi 2 (P)- (32) 
p 

Note that A 2 (£i,£ 2 ) depends on the links £\ and £ 2 but only on their relative orien- 
tation and the distance between them. Of course the eigenvalue also depends on (3 
but we suppress this from our notation. 

Similarly, if the links £\,£ 2 , ■ ■ ■ ,£ n are all distinct, then 

n 

$Q,.. Q „A...a^) = n < w p ) (33) 

1=1 

is an eigenfunction of Kp with the eigenvalue 

n 

A n (£ 1} ...,£ n )=j:e mP) Y[vdP)- (34) 

P i=l 

The eigenvalue is symmetric under permutations of the £j's and simultaneous lattice 
rotations or translations of the links. 

If £ is a link in the lattice T 2 let a\ and a] be the two Ising spin variables on the 
dual sites adjacent to the link £. Then T)e(P) = a\{P)a1{P) where &\(P) is the value 
taken by a\ in the spin configuration corresponding to the loop P. The expectation 
value for an Ising model on T 2 is given by 

((.)) = A 1 E^ <P|0> (-) (35) 
p 
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so 



A n (4,...,4)=Ao(II44>- ( 36 ) 



The correlation inequalities 



n k n 

<IR*i> > <IR4>< II «) (37) 

i=l i=l i=k+l 

now imply 

A n (4,...,4)A > A fc (£ 1 ,...,4)A„„ fc (4+i,...,4). (38) 

We expect that 

A„(A,...,4) < A 1? (39) 

for any n > 1 but do not have a general proof of this inequality It can be checked 
in special cases using the explicit form of the two spin correlation function [^D| and 
the cluster property of correlations. 

4 Multiplicities 

As remarked above, the largest eigenvalue Ao is simple by the Perron-Frobenius 
theorem. Let us consider the first nontrivial eigenvalue A x with eigenfunctions $q/. 
Since is a multiple of by the constant r]i(Q) it suffices to consider the 
functions $o,£ — Vi- We claim that these functions are linearly independent so the 
multiplicity of Ai is at least 2M 2 . 

Suppose there are constants q such that 

£q % (P) = (40) 

e 

for any P e II. We adopt the convention that sums on I run over all links in T2 
unless otherwise specified. Then 

-E<* + E c * = °- ( 41 ) 

By taking P = we also obtain 

5> = 0. (42) 



It follows that 



E ce = (43) 
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for any P e II. Let now x and y be two different lattice points in T2 and suppose 
P is made up of two simple nonintersecting curves 71 and 72 from x to y. Then 
clearly 

£<* = "£<*■ (44) 

£671 ^£72 

Let 73 be one more simple curve from x to y which also avoids 71 and 72. Then 
we can join 73 with either 71 or 72 to make a closed curve (i.e. a loop) in II. We 
conclude that 

£ ci = £ ct = £ c t (45) 

££71 £673 £672 

which implies that 

£ c e = (46) 

for z = l,2, 3. Now take x and y to be nearest neighbours and 71 the curve that 
joins them by one link l\. It follows that C4 = and hence q = for all £ since x 
and ?/ are arbitrary. 

We could have obtained the above result more easily by noting that the functions 
rji are mutually orthogonal in the natural inner product on Ti which we denote by 
( • , • ). This follows from the fact that 774774 with i\ 7^ £2 is an eigenfunction of Kp 
with an eigenvalue A 2 (^i,£2) 7^ A so 77^774 is orthogonal to the constant function, 
i.e. 

(1,77^774) =£77^)774(^ = 0. (47) 
p 

Hence, 

(774,774)= 5 hh 2 M2 (48) 

since the total number of loops in II is equal to 2 M ' 2 . This can also be seen by a 
direct calculation. 

If the inequality (|39|) is valid then the multiplicity of Ai is exactly 2M 2 . This 
follows from the fact that the functions $Q 1 ...Q n j 1 ...e n span Ti. We now prove this 
spanning property. 

Let ^1,^2,^3,^4 be four different links which contain the same vertex. We shall 
call such a collection of links a star. In this case 

r ]il (P)r ]i2 (P)r ]i3 (P)r ]u (P) = l (49) 
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for any loop P since P contains 0, 2 or 4 of the links Ei, . . . ,£4. If £i,£ 2 , £3, £4 are a 
star and n > 4 it follows that 

®Q 1 ...Q n ,l t ..l n = ±®Q 5 ...Q n ,e 5 -tn- ( 50 ) 

Let £ = {£1, . . . , £ n } be a collection of links in T 2 and define 

E c {P) = \{m{P)- (51) 

i&c 

Clearly Ec = ±&Q 1 ...Q n j 1 ...e n so in order to prove the spanning property it suffices to 
identify 2 m2 mutually orthogonal functions of the form Ec- We adopt the convention 
E = l. 

If L and CJ are two collections of links we say that they are equivalent if the 
symmetric difference CAC is a star or the symmetric difference of two or more 
stars. This defines an equivalence relation on the set of all collections of links. If £ 
is a collection of links let [£} denote the equivalence class of C. It is easy to see that 
the total number of elements in each equivalence class is 2 m2 and the number of 
different equivalence classes is also 2 Af2 since the total number of collections of links 
from T 2 is 4 m2 . If we choose one £j from each equivalence class i — 1, . . . , 2 m2 , 
then 

(E Ci ,E Cj )=Yl II ru(P), (52) 
p eedj 

where Cij is a maximal subset of CiACj which contains no star. Hence, for i 7^ j, 

(E Ct ,E Cj ) = (l,E Cij ) = (53) 

since is nonempty for i ^ j. This proves that {E Cl } is a family of mutually 
orthogonal functions. 

5 Discussion 

Using the results derived in the previous section it is quite easy to calculate corre- 
lation function, i.e. functions of the form 

G<P(P,Q)= £ e~^ M \ (54) 

dM=PUQ 

where the loops P and Q lie in two intermediate constant time planes separated by 
N lattice spacings. In Eq. (|54j ) we sum over all surfaces whose intersection with the 
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two intermediate planes are P and Q and S\(M) is a modified action functional, 
counting only edges that are orthogonal to the time direction. We have 

Gf\p,Q) = {8 P ,Kp Q ), (55) 

where Sp and Sq are delta functions in TC. Up to normalization we can interpret 
Gp(P, Q) as the probability of having the loop P at time iV given that we have Q 
at time 0. An easy calculation gives 

2 M 2 

Gf Q) = 2- M2 £ E Ci (PAQ)A», (56) 

where the sum runs over the orthogonal family of functions constructed in the last 
section and A c . is the eigenvalue of Kp corresponding to E^. 

Unfortunately it is not clear how to extend the analysis of the present paper to 
include the curvature terms in the action of the full model. The original transfer 
matrix ([5]) is not invariant under the group Q, only translations and rotations remain 
symmetries. The constant function is not an eigenfunction any more. The only result 
that survives is the fact that the largest eigenvalue is in this case also simple by the 
Perron-Frobenius theorem. 

However, we have managed to solve exactly a three dimensional lattice model 
that describes the diffusion of loops. As a surface model it is not isotropic but can 
be viewed as a stack of two-dimensional Ising systems with an intereaction that is 
sufficiently simple for us to solve the model exactly. 

Acknowledgement. T. J. is indebted to NRCPS "Demokritos" and the CERN 
Theory Division for hospitality. 
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